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Abstract 

Based on the equal-time Wigner function for the Klein-Gordon 
field, we discuss analytically the mechanism of pair creation in a clas- 
sical electromagnetic field including back-reaction. It is shown that 
the equations of motion for the Wigner function can be reduced to 
a variable-frequency oscillator. The pair-creation rate results then 
from a calculation analogous to barrier penetration in nonrelativistic 
quantum mechanics. The Wigner function allows one to utilize this 
treatment for the formulation of an effective transport theory for the 
back-reaction problem with a pair-creation source term including Bose 
enhancement. 



1 Introduction 

The usage of the equal-time Wigner transform of the two-point correlation 
function of a quantum field theory has been proposed to investigate nonper- 
turbatively the time evolution of the vacuum state [I], fj. This approach is 
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based on the Fourier transform of the two-point correlation function split in 
the space coordinates but evaluated at equal times. The equation of motion 
for this object can be formulated for the interaction of the Dirac |jj or the 
Klein-Gordon field 0] with a classical electromagnetic field and exhibits fea- 
tures of a transport equation. It does not, in particular, entail a constraint 
equation — as does the covariant Wigner function || — and can be computed 
as an initial-value problem. This makes this approach a good candidate for 
problems like the formation of the quark-gluon plasma and the back reaction 
from pair creation in the early universe. 

The problem of back-reaction from pair creation in a classical electromag- 
netic field has been investigated earlier using the field equations directly [§]. 
Numerical results ]5|, |(J indicate that the exact solution of the field equations 
can be very well approximated by a transport theory involving a Schwinger 
term and Bose enhancement or Pauli blocking, as pertinent. This has raised 
the question of how to derive such a theory consistently starting from the 
field equations. While this has been attempted before using the covariant 
Wigner function , we here show how the pair-creation term can be isolated 
beginning from the equation of motion for the equal-time Wigner function of 
a Klein-Gordon field interacting with a classical electromagnetic field. Simi- 
lar approaches based on the equal-time Wigner operator have been taken in 



2 Formalism 



2.1 The equal-time Wigner function 

The Wigner function of the Klein- Gordon field can be defined from the sym- 
metrized two-point function 

C+ p (q,q';t) = (n\{$ a (q,t),$+(q',t)}\n) (1) 

evaluated at equal times t (q and q' are three-dimensional vectors). In or- 
der to obtain an evolution equation for the Wigner function, it is necessary 
to consider the field as a two-component object $ in the Feshbach-Villars 
representation 

*=(tl> = ^ + X, (i^--eA°) ( J ) = m^-x), (2) 
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where designates the basic Klein-Gordon field obeying the equation of 
motion in a classical electromagnetic field A^, 

- ieAfj) (0" - ieA") + m 2 ) <f>{x) = 0. (3) 
The Wigner function is defined by 



f ( % ie fi+y/z \ 

P(q,P) = J d3 2/ C + (q- \y,q+\y,t) expi-p-y + — J A{x) ■ dx J , 

(4) 

where A(x) is the electromagnetic vector potential. The equation of motion 
of the Wigner function follows from the equation of motion for the Feshbach- 
Villars field 

with the Hamiltonian 

H= ^ ~ eA)2 a + mb + eA°l . (6) 
2m 

Here p designates the derivative operator —id/dq, and a and b are the ma- 
trices 

a = (<7 3 + iV 2 ) = \ , & = ^3 = J — 1 ) ' ^ 

In the free case, (0) has the solutions 

$(q,t) =u i±) (p)e Tl{E - t - p - q \ (8) 
with 

representing particles of positive and negative energy, respectively. 

Applying phase-space calculus 0, the equation of motion for P can be 
written as a differential equation in phase space, 

( d d \ 
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p ® _|_ e (^S(g) x ^ ) ^ + 
2m <9g V mj dp 



(a • P(g,p) + P{q,p) ■ a + ) 




1 d 2 \ 

P(q,p)-P(q,p)-b), (10) 

where the dots indicate quantum corrections involving higher derivatives of 
P and of the electric and magnetic field. 

The first term on the right-hand side is a (nonrelativistic) flow term. The 
second and third terms incorporate relativistic effects (manifested by the 
appearance of p 2 and m 2 ) and interferences on the Compton scale which are 
present in any Wigner function and result from the impossibility to measure 
both position and momentum of a quantum-mechanical particle at the same 
time. Since the equation of motion for the Wigner function is exact, it 
incorporates a mechanism to generate these interferences — namely the term 
d 2 /dq 2 — which can be safely neglected in the semiclassical limit, i.e., when 
the spatial variation of the Wigner function is small on the Compton scale. 
Note that we do not have to make any assumption about the variation in 
time of the Wigner function. This feature is special to the equal-time Wigner 
function. This can also be seen from the fact that, after expanding P(q,p) on 
a matrix basis (see [[| and section |2.2j ), p 2 appears only in the combination 

1 d 2 

p 2 + m 2 - -— ^p 2 + m 2 = El. (11) 
4 aq l p 

With this approximation, which is exact for a spatially homogeneous prob- 
lem, eq. (pTOf ) contains both flow terms as well as local couplings involving m 
and p 2 . In particular, the equations have nonstationary solutions even in the 
case of a spatially homogeneous Wigner function. These oscillations are due 
to the internal charge degree of freedom of the Wigner function. 



2.2 Expansion in a diagonalizing basis 

To analyze these local oscillations, we expand P into eigenoscillations. In 
the free case, two of the eigenoscillations have the frequencies ±2E p , while 
two others have zero eigenfrequency and correspond to constant solutions 
of the equations of motion. The latter can be interpreted as semiclassical 
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quantities — which later may serve to define the physical phase-space densities 
of particles and antiparticles — while the former represent interferences. 
Expanding into eigenvectors of the local oscillations, 



i=i 



P(q,p) = J2fi P n (12) 

with the p-dependent basis 



1 / m/Ep + Ep/m m/Ep — Ep/m 
4 I m/ E p — Ep/m m/E p + E p /m 



P3/1 



^(uWgnW + uH^nH), (13) 

U I J ) = I ® ^ (+) - «H ® «<">) , (14) 

if m/Ep-Ep/m m/E p + E p /m =p 1/2 \ 
4 \ m/^ - Ep/m ±1/2 m/Ep - E p /m J 

M (T) ®M (±) , (15) 



the corresponding semiclassical equations of motion for / are 



(16) 



E 2 2" p ^ 3 ' ( 8 ^ 

(l +e ^ } -|) /4 = ^l;-| /2 

S(q)-p fi 

E 2 ~2 ~ p ^ 4 ' ^ 9 ^ 

Here the approximation of eq. (|ll|) has been made, and the higher-order 
quantum corrections to the electromagnetic interaction have been neglected. 
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The flow term now has the correct relativistic form. The self-coupling terms 
that give rise to local oscillations are isolated in the third and fourth equation 
(which are complex conjugate to each other), rendering f 3 and / 4 oscillatory- 
even in the case of a homogeneous Wigner function. The expansion into 
shows that the nonoscillating terms do not mix positive- and negative-energy 
solutions, while the oscillating ones do. The latter oscillations are therefore 
caused by Zitterbewegung interferences between particles and antiparticles 
and would not be present in a nonrelativistic formulation. 

The phase-space densities of observables can be expressed in these vari- 
ables: 

charge: ef 2 , 
energy: E p f lt 

current: (fx + f 3 + / 4 ) , 

momentum: — pf 2 + ip{fz — ft)- (20) 

Note that the current involves the highly-oscillatory components f 3 and / 4 . 
It couples to the electromagnetic field and thus provides a mechanism for the 
emission of photons in particle annihilation. 
A new self-coupling term with the coefficient 

eJi^P (21) 

has appeared. This term is due to the fact that the basis fll3|) depends on the 
coordinate p and therefore generates an additional term when acted upon by 
a force term. The physical reason for this term is that the local oscillations 
fs and / 4 have different "directions" in the /-space at different momenta of 
the particle. When a particle changes its momentum because of an external 
force, its local oscillations will be slightly out of phase and arouse the f\- 
component, which can be interpreted as pair creation. Inspection of (|13|) 
reveals that the basis changes most at low momenta, and thus the coupling 
term is maximal at p <C E. Physically this means that the pair creation 
happens at approximately zero momentum of the particles. 

The equal-time Wigner function is, by definition, not Lorentz-invariant. 
In another reference frame, the Wigner function involves the expectation val- 
ues of the fields at different times. It is not clear that the new function can 
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be constructed from the Wigner function in the original frame of reference 
since phase information may be lost in the equal-time Wigner transforma- 
tion. Moreover, since our description is valid in any inertial system, it leads 
to the apparent paradox that pair creation happens at approximately zero 
momentum in any frame of reference. A similar paradox is discussed in ||. 

The understanding of these equations can be further facilitated by con- 
sidering the linear combinations 

/+=|(/l + / 2 ), /i = /+ + /_, 

/- = !(/i -/a), h = /+-/-, (22) 

which give the phase-space densities of particles /+ and antiparticles /_. 
Rewriting the above equations yields 



K dt ^ J dp E p dq) J ^ 2E p dq 

+ e^(fs + U), 

'Z + *l*)-l-4r-l)f- - ^ItA + A) 



dt dp E p dq J 2E P dq 



£ m p 

- e ^^(/s + fa), 

d d \ 1 p d 

Jt +e£{q) -dp-) h = -2E v -dq h 



El 2 

[i +e£{q) -lp) u = 4i-| /2 

+ eS-^-2iE p f 4 , (23) 

which exhibits clearly a flow term for / + and /_ on the left-hand side. The 
components f% and f'4 do not have a flow term for position space since the 
relevant interferences are not localized. The right-hand side contains quan- 
tum effects and, in particular, pair creation. In the nonrelativistic limit, /+ 
and /_ change on a time scale much larger than the frequency of f% and 
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fi-, so that the right-hand sides average out and the equations describe a 
relativistic free gas of positively and negatively charged particles. Note that 
f-(p) refers to antiparticles of momentum — p which is consistent with the 
usual hole interpretation. 



3 Pair creation 



3.1 Reduction to a variable- frequency oscillator 

Let us now investigate the case of a homogeneous constant field, that is, 
the Schwinger mechanism of spontaneous pair creation. Because of the spa- 
tial homogeneity, the approximations on local interferences made above in 
eq. fllQ|) hold exactly. We shall thus be able to derive exact equations for the 
Wigner function in this case. 

We first remove the force terms from the equations of motions by sepa- 
rating out the classical motion using the method of characteristics. Consider 
a new function / (p , t) defined by 

l(p,t) = J dp l(p ,t)8(p-p(p ,t)), (24) 

where the function p(po,t) describes the classical motion of a particle with 
initial momentum po, 

p{po,t)= Po + [ dt'e£{t). (25) 

J — oo 

The phase-space evolution of / (p, t) translates into 



(! +e£ '£U< p - ( >=/ dp » 



%^( P -p(p„, ( )). (26) 



In this way, the momentum derivative is absorbed into the classical evolution. 
Here f (po,t) represents the fate of a phase-space element that started at 
p = Po and moves in phase space according to Newton's law. It may therefore 
be called a test particle whose inner degrees of freedom in the Feshbach- 
Villars sense are governed by 
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dt 



0, (28) 



dt E 2 p 2 

dh eS-p h 
dt El 2 



+ 2iE p f 3 , (29) 
- 2iE P h, (30) 



where p = p(po, t) is now time-dependent. The set of partial differential equa- 
tions (p~3| ) is thus reduced to a set of ordinary differential equations. Instead 
of solving the flow equation in phase space directly, we have in this way intro- 
duced a 'test particle', moving in phase space according to Newton's law and 
possessing an internal degree of freedom governed by (f27[)-(|3"0"l). Accordingly, 
the equations must be solved for all possible initial values Pq. The function 
/ (po, t) thus computed can then be assembled using (j24j) to give the actual 
solution of the system of partial differential equations (|TB|). 

Due to spatial homogeneity, the charge density vanishes. Initially, only 
fx will be nonzero, but by the action of the electric field, the interference 
components f 3 and / 4 are aroused which in turn couple back to fx. The 
coupling 

e m ( 31 ) 

V 

vanishes as p —>■ oo. At large times, fz and / 4 will therefore settle down into 
steady oscillations of increasing frequency and decouple from fx. 

The set of ordinary differential equations (E7p is equivalent to a single 
differential equation in a new variable £(t), 

n2 



£ + = 0, 



through the substitution 



fi 

h+h 



E p 



24 0) V ' El 



2Ei 0) V ' El 



1 





(32) 


lei 2 ), 


(33) 


\i\ 2 ) , 


(34) 




(35) 
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where dots denote derivatives with respect to t, and E^ is the test-particle 
energy at some fixed reference time. The relationship between pair creation 
in a homogeneous field and oscillators of time-dependent frequency has been 



noticed much earlier from the wave-function approach to pair creation [10 



11, |T|, 13] . By using the solutions given in [12|], we will be able to give an 



approximate asymptotic solution for the Wigner function. 
3.2 Approximate solution 

The general solution of the variable frequency oscillator (RI3) has the form 



m = ~^=C 1 e- ir + -^=C 2 e iT , (36) 

where G\ and C*2 are slowly-varying functions, and r(t) is defined by 

r(t) = f&t'E p {t'), (37) 

with an arbitrarily chosen reference point t\. Making the Liouville substitu- 
tion 

m = 4=x(r(t)), (38) 

one obtains the oscillator equation in x( T ) 

X " + (l + h(r))x = 0, (39) 
where primes denote derivatives with respect to r, and the function h(r) is 

Equation (^) is the time-dependent Schrodinger equation of a particle in a 
potential h{r). The shape of the potential is determined by E p (r) and thus 
by the external forces. In particular, h(r) vanishes for times r when E p (t) 
is constant, i.e., when the external field is switched off. Even in a constant 
electric field, E p ~ r for t — ► ±oo, r ~ t 2 and thus h(r) ~ r -2 . The potential 
is therefore localized, and the problem is equivalent to a barrier-penetration 
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problem in quantum mechanics. Note that for a constant electric field S, the 
barrier potential is 

Mr) 4 fe^ 2 . (41) 



4\ f? v 



in which the quantity (pT[) resurfaces. 

The coefficients G\ and C 2 represent the incoming and the reflected waves, 
respectively. They are asymptotically constant, and the process is character- 
ized by a barrier penetration coefficient (in addition to a phase shift) 



Po 



Co 



Ci 



2 



t -> ±00, (42) 



the index indicating a transition from the vacuum value. Since this prob- 
lem is similar to the problem of barrier penetration in ordinary quantum 
mechanics, the calculation of this parameter is accomplished by well-known 
methods. The quantity po will depend on the initial momentum and on the 
electric field. 

Substituting eq. (06|) and r = E p t into the relations ( [33| - |35D , we obtain 

A = iCil 2 + IC 2 [ 2 = |Ci[ 2 (l + po), (43) 
f 3 = 2C*C 2 e 2iEpt , (44) 
U = 2C*C ie - 2iEpt . (45) 

For the vacuum, we choose 

d = 1, C 2 = 0. (46) 

The relationship between the quantity p and the number of particles can be 
obtained by considering that, due to the reality of the oscillator equation, 
the form — ££* is a constant, 

Ci ~ iC = -2iE p (\d\ 2 + \C 2 \ 2 ) = const = -2iEf\ (47) 

where at t —>■ — 00, C\ = 1 and C 2 = have been chosen. From p = \C 2 /Ci\ 2 , 
it now follows that 

1 1 EM 

Idl 2 = 1 \C 2 \ 2 = p p „ . (48) 

' ' l-pE p (tY 1 21 H l-pE p (t) 1 ; 
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In particular, we get 



£ (0) 1 + o 

^ 2 + ^ 2 -W)^ m 



which is related to the resulting current density. Substituting fl3q ) into 
it follows that 

h = fHiAi 2 +ic 2 i 2 )=^, (50) 

h + U = -4y(qrC 2 e-*+c.c), (51) 
2 hip 

i(h-U) = -%y (l^l 2 + |C 2 | 2 - CrC 2 e- 2i - - c. c.) . (52) 
AEp ' v ' 

As was to be expected, j\ is asymptotically constant, while the other two 
components oscillate because of interferences. The current density in phase 
space is 

3(q,p) = ^-(fi+fs + h)- (53) 

The components and J4 will give rise to fast oscillations. Since the current 
couples to the electromagnetic field, this may cause the emission of photons 
(or, in the semiclassical approach, electromagnetic waves of frequency 2E p ). 
For our purposes with coupling to a classical electromagnetic field, the oscil- 
lations tend to cancel out upon measurement, and the main contribution to 
the current density in phase space comes from 

The vacuum value of this quantity is ep/E p ; it represents the charge of the 
Dirac sea and must be subtracted to yield a measurable quantity. Thus the 
phase-space density of pairs produced is given by 

n=j^-, (55) 
1 - p 



in full correspondence with eq. (3.7) of [12 
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This yields the following picture of pair creation: Let us assume a constant 
field £ switched on at some time to. For each initial momentum po, we have 
at t — > — oo G\ = 0, C 2 = 0, and thus fx = 1, / 3 = / 4 = 0. When 
the field is switched on, the particle is accelerated or decelerated. In the 
latter case it reaches its minimum velocity at some time, corresponding to 
the barrier in the equivalent quantum mechanical problem. In this moment, 



C\ and C 2 change so that fx, f 3 , and / 4 acquire the values (p0|) - (j52|) . This 
transition takes place when p(po,t) is minimal, i.e., when its parallel part 
p\\ vanishes. Looking at the phase-space distribution / (p, t) we therefore 
see the characteristic trough found in M and which signals the impending 
pair production at a constant rate per unit space and time: The phase-space 
density changes at p\\ = which creates the left end of the trough; the change 
is then transported away by the electric field, and the edge that appeared 
when the field was switched on constitutes the right end of the trough. 

When the field is not constant in time, it may happen that a test particle 
reaches p\\ = more than once. We then can observe the effects of Bose 
enhancement (or Pauli blocking for fermions). In a very strong electric field 
with appreciable pair creation we can assume that the transition takes place 
during a very small time. This is comparable to solving the barrier pene- 
tration problem by gluing together plane-wave segments at the barrier with 
different amplitudes and phase shifts. The time evolution of £ (t) can then be 
thought of as being described by ( j36| ) with the coefficients C\{t) and C^i) 
changing their value abruptly whenever p(t) is small. Let us assume that any 
one of these transitions can be described by 

e -iE p t _^ A e -i(E p t+a) + B e i{E v t+f5) ^ ( 5g ) 

where the amplitudes A, B and the phases a, j3 depend on the field strength 
at the transition. If a general initial state is given by 

with some positive-state admixture characterized by p and a phase shift S, 
this initial state will go into a new state 

(Ae~ ia + ^e^) e~ iE ^ 
+ (Be^ + ^7le ,(a+i) )e^. (58) 
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The new positive-state admixture p', defined as the ratio of the square of the 
amplitudes in front of the oscillatory terms, is accordingly given by 



p + p + 2 v /pp^cos(a + 5 - (3) 

1 + PoP + 2 v /pp^cos(o; + 6 - p) 

where p = B 2 /A 2 is a characteristic quantity of the transition and will 
depend on the electric field strength at the time of the transition. In a 
nearly constant field, p will be given by the Schwinger term. As we see, 
p behaves as an additive quantity in the transitions, up to an interference 
term. The latter can be assumed to average out when different phase-space 
trajectories are considered. 

Because of this, we can assume that the time evolution of the quantity 
p(t) is given by 

p(t)= f dt'6(p(p ,t'))p (t'), (60) 

where p(po, f) is the momentum of the test particle (pit) depends of course on 
the parameter p ), and po(t') is determined by the electric field at the time t'. 
As a result of this assumption, p(t) will be a sum of all contributions coming 
from the individual transitions which occur whenever p(po,t') is small. The 
approximation here lies in neglecting the actual transition process in which 
p changes to p + p , since the duration of the transition can be assumed to 
be quite small in strong fields. 

3.3 Constructing a flow equation with source term 

We have thus derived an expression for the phase-space density of created 
pairs. Remembering that p = p(p ,t) we now consider the phase-space func- 
tion 

n( P0 ,t) = / {P f tV ^ 
1 - P(Po,t) 

where the tilde reminds us that the classical motion of the test particle has 
been separated out. To transform back to a true phase-space equation, we 
take the time-derivative of n(p ,t), 

dn(po,t) 1 dp 2 dp 

^^ = JT^pW^ = {1 + nM) W (62) 
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\eS(t)\8{p(p ,t))p (t). (63) 



Substituting the time evolution of p(t), eq. fl6"0|), this is 

dn(po,t) 
dt 

Defining the true phase-space function n(p,t) as in (Hj), 

n (p, t) = J dp n(p , t)5(p- p(po, £)) , (64) 
the phase-space equation of motion is 

(j t +e£-^jn{p,t) = (l + n(p,t)f\e£(t)\ 

x / dpoPo(t)5(p(p ,t) 5(p-p(p ,t)) 

= (l + n(p,t)) 2 Po (t)5(p). (65) 

This is a phase-space flow equation with a Bose-enhanced production term on 
the right-hand side. The quantity po(t), characterizing the pair-production 
rate, can be computed either perturbatively (pair production by oscillat- 
ing fields) or nonperturbatively (Schwinger term), or by some combination, 
yielding an enhanced estimate of the pair-production rate. 

It has been found earlier || in the study of the back-reaction problem 
that pair creation in a spatially homogeneous problem can be surprisingly 
well described by the source term 



;i + 2n(p,t))ipin 



1 + exp 



Tim 2 



S(p) 



I €-S I / TCTfi^ \ 

w (l + 2n(p,t)y-^expl- — j 5{p) , (66) 

the difference in the Bose enhancement factors in Eqs. (|6"5|) and ( |56] ) stemming 
from the lack of a mechanism in the former that annihilates pairs in favor of 
an augmentation of the electric field. 



4 Summary 

We started from the equation of motion for the Wigner function. This equa- 
tion of motion is a direct consequence of the field equations. We have then 
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derived an approximate solution in terms of a transport equation involving a 
Schwinger pair creation and a Bose enhancement term. The Wigner function 
exhibits an internal structure causing local oscillations, i.e., nonstationary so- 
lutions in a spatially homogeneous situation. The form of these oscillations, 
corresponding to the internal charge degree of freedom of the Klein-Gordon 
field, are different at different momenta of the particles. Thus in the pres- 
ence of an external force the oscillations mix with the physically observable 
phase-space densities of particles and antiparticles and give rise to particle 
creation. This process is described analogously to a barrier-penetration prob- 
lem in quantum mechanics. The semiclassical approximation enables us to 
calculate the pair creation rate and to give an estimate for repeated pair 
creation involving Bose enhancement. 

While it was known earlier that pair creation in an external field can 
be described by a variable-frequency oscillator from considering the wave 
functions, it is the advantage of the Wigner function that it allows one to 
cast this in a phase-space formulation. The use of the equal-time Wigner 
function is here especially advantageous over the split-time Wigner function 
since it allows us to exploit the time- dependent problem of the variable- 
frequency oscillator and thus the known solutions to it in a very natural 
way. 

Some open questions remain: For spatially inhomogeneous fields, we have 
only shown how the homogeneous-field limit can be recovered. By making 
a systematic expansion in the field gradient, it should be possible to derive 
corrections to this result. We have also restricted ourselves to classical elec- 
tromagnetic fields. In a next step, the electromagnetic field should be treated 
dynamically. In this approximation, emission of photons, pair annihilation 
into photons, and particle interaction enter. Finally, even in the semiclassical 
limit, the theory is, in general, in need of renormalization 

In its current form, the equal-time Wigner function is especially suitable 
for the explicit solution of semiclassical problems like particle creation in def- 
inite boundary conditions. The fact that it is not relativistically invariant 
makes it more difficult to derive generic features of quantum field theories 
from it. Still, a possible extension could be to consider general iV-particle 
Wigner functions constructed from 2N field operators. The resulting hierar- 
chy of equations can then be systematically studied. 
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A Appendix: The Schwinger term, obtained 
by means of the imaginary time method 



Since fl3"2| ) corresponds to the time-independent Schrodinger equation de- 
scribing tunneling through a potential barrier, p is the reflection coefficient 
for this problem and can be computed by a standard method, namely that 
of imaginary times |T2, I3fl . In this case the solution is determined by the 



singularities of h(r) in the complex plane: If h(r) has its singularity at To, 
the value of the barrier reflection coefficient is 

p = 4cos 2 (7r/2a) exp(-4Imr ), (67) 

where a is deduced from the behavior of E p (t) near the singularity, 

u(t) ^LU (t-t ) a -\ (68) 

In the case of a constant field, 



E p (t) = ^(p + eSt) 2 + m 2 } (69) 
and the singularities are at E p = 0, corresponding to 

t0 = za±p, (70) 

yielding a = 3/2. To compute the corresponding value of r , we choose 



the integration contour in (|37D along the real axis from t\ to Re to and then 
parallel to the imaginary axis to Imt - No contribution to Imr arises along 
the real axis, so the integral to do remains 



/■Imio 

Imr = Im i / d(Imt) EJt) 
Jo 



(71) 

Re t=Re to 
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Writing x = Im£ , this is 



rim to 

Imr = Re / dxielS 
Jo 



x 



x 2 + 2i[ Re t + 



m , 



P\\ 



e\£\l e 2 |^l 2 



1/2 



= T7T 



ml 



4e|£| 



The barrier-reflection coefficient is therefore 



(72) 



m i 



which agrees with the lowest-order term of the Schwinger formula. 



(73) 
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